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ABSTRACT 


An  experimental  investigation  was  performed  to  evalu¬ 
ate  tearing  instability  theory  by  varying  the  applied 
tearing  modulus,  Tappiiecj,  so  that  fracture  instability 
would  be  initiated  at  various  levels  of  stable  crack  exten¬ 
sion.  This  is  an  extension  of  past  investigations  of 
tearing  instability  theory  in  that  crack  extension  was 
monitored  continuously  using  the  D.C.  potential  drop 
technique,  enabling  the  applied  and  material  tearing  moduli 
to  be  calculated  at  the  point  of  instability. 

The  results  of  this  investigation  indicate  that,  in 
most  cases,  fracture  instability  occurred  when  the  difference 
between  the  applied  and  material  tearing  moduli  was  on  the 
order  of  10Z.  Variations  in  the  load  versus  displacement 
records  of  the  specimens  near  maximum  load  due  to  local  in¬ 
stabilities  and  friction  in  the  load  train  precluded  measure¬ 
ment  of  a  smooth  applied  tearing  modulus  curve. 
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INTRODUCTION 


The  tearing  instability  concept  was  developed  by  Parish  as  a  means  of 
assessing  the  stability  of  flawed  structures  that  fail  by  ductile  tearing. 
Tearing  instability  theory  states  that  a  flawed  structure  loaded  beyond  Jjc 
and  at  limit  load  will  tear  in  a  stable  manner  so  long  as  the  applied  tearing 
modulus,  TappHe(j,  is  less  than  the  material  tearing  modulus,  Tmaterfai. 

Tearing  will  proceed  in  an  unstable  manner  whenever  the  ^applied  equals  or 
exceeds  Traaterja^.  Several  experimental  investigations  were  conducted  to 
validate  tearing  instability  theory  for  materia' s  exhibiting  both  high  and  low 
toughness  [2-4].  Vassilaros,  et  al.^  used  the  compact  tension  specimen  in  a 
variably  compliant  screw-driven  test  machine  to  validate  the  tearing  instability 
concept  for  several  high  toughness  alloys  with  Tmater^a^  spanning  the  range  of 
12  to  170.  The  Tapp^iecj  formulations  used  to  analyze  the  test  data  included  the 
generalized  expression  of  Paris,  et  al.l,  a  modified  Paris  formulation,  and 
an  expression  developed  by  Ernst,  et  al.^  using  the  key  curve  approach.  Among 
the  conclusions  of  the  investigation  were  that  the  Tappued  expression  of 
Paris,  et  al.  was  non-conservative  in  predicting  the  point  of  instability  in 
the  compact  specimen  and  that  the  generalized  Paris  analysis  evaluated  at 
maximum  load  was  most  consistent  in  predicting  instability.  They  also  concluded 
that  the  analysis  of  Ernst,  et  al.  appeared  to  be  accurate,  but  the  analysis 
required  precision  beyond  that  which  was  available  in  their  test  data. 

The  purpose  of  the  present  investigation  was  an  extension  of  the 
Vassilaros,  et  al.^  study  to  further  evaluate  tearing  instability  theory  by 
varying  Tappne(j  so  that  instability  would  be  initiated  at  various  levels  of 
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stable  crack  extension  and  to  employ  the  Ernst  analysis^  to  calculate 
Tapplied • 

EXPERIMENTAL  PROCEDURE 
Material 

The  material  used  for  all  tests  conducted  in  this  investigation  was 
ASTM  A106  steel.  The  specimens  were  removed  from  a  39  inch  (990  mm)  diameter 
pipe  with  a  3.75  inch  (95  mm)  thick  wall.  The  specimens  were  tested  in  the  L- 
C  orientation.  The  chemical  composition  of  this  steel  is  presented  in  Table  1 
and  typical  room  temperature  mechanical  properties  are  listed  in  Table  2.  A 
Charpy  energy  -  temperature  transion  curve  produced  from  specimens  oriented  in 
the  L-C  direction  is  plotted  in  Figure  1.  All  J-integral  resistance  curve 
tests  were  conducted  at  room  temperature  which  is  in  the  transition  region  for 
this  steel. 

Testing  Procedure 

The  specimens  were  1  inch  (25  mm)  thick  IT  compact  tension  specimens 
with  20%  side  grooves  as  shown  in  Figure  2.  The  specimens  were  fatigue 
precracked  to  a  nominal  initial  crack  length  of  aD/W  =  0.65.  The  tests  were 
conducted  in  a  screw-driven  test  machine  with  a  spring  in  series  which  allowed 
the  compliance  of  the  testing  machine  to  be  varied  between  tests,  Figure  3. 

The  spring  consisted  of  a  pair  of  1  inch  (25  mm)  thick  titanium  plates 
separated  by  a  pair  of  rollers.  By  changing  the  span  of  the  rollers,  a  range 
of  stiffnesses  could  be  achieved;  a  short  span  produced  a  stiff  system  and  a 
large  span  produced  a  relatively  compliant  loading  system.  This  is  the  same 
system  used  by  Joyce  and  Vassllaros^,  and  Vassilaros,  et  al.^  in  their  previous 
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investigations.  The  stiffness  of  the  testing  machine  was  varied  between 
10500  lbs/in  (1840  N/mm)  and  7400  lbs/in  (1300  N/mm). 


During  the  tearing  instability  tests,  load,  crack  opening  displacement, 
and  system  compliance  were  recorded  digitally  at  approximately  0.5  s  intervals 
with  an  interactive  computer  system.  The  crack  length  was  monitored  contin¬ 
uously  by  using  the  D.C.  potential  drop  technique.  High  speed  digital  oscill¬ 
oscopes  were  also  used  to  record  the  load,  crack  opening  displacement,  compliance 
and  crack  length  during  the  instability  event.  The  data  acquisition  rate  during 
the  instability  event  was  1  ps  per  point.  The  oscilloscopes  were  triggered 
by  the  sudden  increase  in  crack  opening  displacement  which  occurred  when  the 
crack  became  unstable.  Pre-triggering  was  used  so  that  data  prior  to  the 
instability  could  be  matched  up  with  the  data  taken  interactively  with  the 
computer . 

ANALYSIS  PROCEDURE 

J  integral  resistance  curves  were  calculated  for  each  test  according  to 
the  expression^ 

JI  =  [Ji  +  HAi ,i+l/biBN] [ l-Y(ai+i-ai)/bi]  (1) 

where 

n  =  2  +  (0.522)bi/W  for  the  compact  specimen 

W  =  Specimen  width 

Y  =l+(0.76)bj/W 

b[  =  instantaneous  length  of  the  remaining  ligament,  (W-a^) 

a|  =  instantaneous  crack  length 

B^  =  minimum  specimen  thickness 

A( t i+2  =  Area  under  the  load  versus  load-line  displacement  record 

between  lines  of  constant  displacement  at  points  i  and  i+1 

The  crack  length  was  determined  from  the  D.C.  potential  drop  data  and 

calculated  according  to  the  expression^ 

a  =  [(PD  -  0.605)0*345/1, 5]  (2) 
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where 

PD  =  (PD-P1)/(P2-Pi) 

Pl,P2  =  calibration  constants 

Equation  (2)  requires  a  two-point  calibration  which  relates  the  initial  and 
final  crack  lengths  to  the  repective  potentials.  For  several  of  the  tests, 
the  instability  event  led  to  complete  separation  of  the  specimen  halves  and  it 
was  not  possible  to  determine  the  calibration  constants.  Pi  and  P2  using  this 
approach.  For  these  cases,  a  value  of  Pi  =  0.2  was  assumed.  This  is  an 
average  value  of  P^  calculated  from  many  tests  conducted  on  IT  compact  specimens 
from  the  same  section  of  A106  steel  pipe.  It  should  be  noted  that  this  value 
of  P^  is  specific  to  the  instrumentation  and  procedure  employed  in  this  testing 
program  and  does  not  apply  in  general  to  other  materials  or  experimental 
apparatus . 

The  material  tearing  modulus  was  calculated  from  the  expression, 

^material  =  (E/oG2)(dJ/da)  (3) 

where 

E  =  elastic  modulus 

a0=  flow  stress  (ay  +  au)/2 

dJ/da  =  slope  of  tne  J-R  curve 

dJ/da  was  calculated  from  a  least-squares  power  law  fit  to  the  J-R  curve  of 
the  form^ 

A2 

J  =  Ai(Aa)  (4) 

The  applied  tearing  modulus  was  calculated  using  the  expression  of  Ernst, 


Y  =  1  +  (0.76)b/W 

b  =  remaining  ligament,  (W-a) 

n  =  2  +  (0.522)b/W 

=  machine  stiffness 
W  =  specimen  width 

E  =  elastic  modulus 

oQ=  flow  stress,  (ay  +  ou)/2 
H'/H  =  (1/P)  (dP/d6)  +  ( n2P/b2 ) [ (dJ/da)  +  yj/b] 
dP/d<5  =  slope  of  the  load-load  point  displacement  curve 
dJ/da  =  slope  of  the  material  J-R  curve 

In  order  to  reduce  oscillations  in  the  calculation  of  Tapplied*  the  load- 

displacement  data  was  smoothed  by  averaging  the  data  over  30  point  intervals. 

The  smoothing  had  a  negligible  effect  on  the  calculation  of  J,  but  provided 

for  a  great  reduction  in  the  oscillations  in  TappHecj  by  providing  smoother 

derivatives  in  the  H^/H  function. 

RESULTS  AND  DISCUSSION 

A  summary  of  the  results  for  the  tests  is  presented  in  Table  3.  The 
power  law  fits  to  the  J-R  curves  from  these  tests  are  plotted  in  Figure  4. 
There  is  a  considerable  amount  of  variability  in  the  J-R  curves  from  specimen 
to  specimen  that  is  not  normally  expected  since  all  of  the  curves  were 
measured  using  the  same  material,  specimen  geometry  and  nominal  initial  crack 
length.  This  variablity  may  be  due  in  part  to  the  location  of  the  specimens 
in  the  wall  of  the  pipe  or  the  fact  that  this  material  is  in  transition  at  the 
test  temperature.  The  difference  in  the  J-R  curves  beyond  Jjc  made  it 
difficult  to  separate  the  effect  of  varying  the  amount  of  stable  crack  growth 
by  changing  the  system  compliance  since  the  material  resistance  to  tearing 
from  specimen  to  specimen  was  not  consistent. 


Table  1.  -  Chemical  Composition  of  ASTM  A106  Steel 


C 

S 

P 

Si 

Ni 

Mn 

Cr 

Fe 

vt  Z 

0.26 

0.017 

0.028 

0.022 

0.26 

a 

0.11 

bal 

Table  2.  -  Mechanical  Properties  of  ASTM  A106  Steel  at 

Room  Temperature 
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Flow 

2 

Reduct . 

J 

Strength 

Strength 

Stress 
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in  Area 
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50.0 

79.0 

64.5 

31  1 

70 
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p  Table  3.  -  Summary  of  Fracture  Test  Results 


Specimen 

ID 

a/W 

Machine 

Stiffness 

lbs/in 

(N/mm) 

J  at 
Inst . 
in-lb/ in2 
(kJ/m2) 

Delta  a 
at  Inst, 
in . 

(mm) 

^material 

^applied 

GAP- 2 

0.68 

8650 

3062 

0.107 

76 

71 

(1514) 

(536) 

(2.72) 

GAF-3* 

0.65 

7850 

3068 

0.077 

125 

123 

(1374) 

(536) 

(1.96) 

GAF-4 

0.65 

9750 

3349 

0.81 

151 

110 

(1710) 

(586) 

(2.04) 

GAF-5* 

0.65 

7400 

3363 

0.094 

120 

85 

(1296) 

(590) 

(2.39) 

GAF-6 

0.65 

10350 

3178 

0.100 

128 

127 

(1812) 

(550) 

(2.54) 

GAF-12 

0.66 

10350 

2813 

0.112 

91 

84 

(1812) 

(490) 

(2.84) 

Fracture  mode  converted  to  cleavage  after  ductile  instability 


A  typical  load-displacement  record  from  specimen  GAF-6  is  shown  in 
Figure  5.  The  data  beyond  the  point  of  instability  was  recorded  on  the 
digital  oscilloscope.  The  value  of  Tmateriai  at  the  point  of  instability  was 
107  and  Tapp^^e£j  was  106,  Figure  6.  These  values  are  in  good  agreement  at  the 
point  of  instability,  but  this  was  not  the  case  for  all  of  the  tests.  An 
example  of  a  test  in  which  instability  occurred  well  before  T^t-gr^i  was 
equal  to  Tapp^^ecj  is  specimen  GAF-4,  Figure  7.  In  this  case,  Tmaterial  was 
126  at  the  point  of  instability  and  Tapp^^ecj  was  only  92,  which  would  have 
predicted  a  stable  response.  In  most  cases,  instability  occurred  when  the 
difference  between  Tapp^ie<j  and  Tmater^aX  was  on  the  order  of  10%,  although 
Tmaterial  was  alwaYs  greater  than  Tappue^  at  instability. 

Several  comments  are  in  order  regarding  the  calculation  of  Tappiied  using 
the  analysis  of  Ernst,  et  al.  As  can  be  seen  from  Figures  6  and  7,  the 
^applied  curve  is  not  a  smooth  function  as  would  be  expected.  This  is  due  to 
the  calculation  of  the  slope  of  the  load-displacement  curve  from  experimental 
data.  All  of  the  functions  that  go  into  the  calculation  of  Tapp are 
smooth  functions  of  the  crack  length  or  crack  extension  with  the  exception  of 
the  slope  of  the  load-displacement  curve  which  is  determined  from  the 
experimental  data.  Local  instabilities  such  as  the  growth  and  coalescence  of 
voids  by  "fast  shear”  and  other  micro-mechanisms  cause  small  fluctuations  in 
the  load  which  lead  to  large  changes  in  the  calculation  of  the  slope  of  the 
load-displacement  record  when  the  curve  exhibits  a  "flat"  nature  as  was 
observed  in  all  tests  conducted.  In  the  vicinity  of  maximum  load,  the 
calculation  of  the  slope  is  strongly  influenced  by  the  resolution  of  the 


instrumentation  and  the  noise  level  present  in  the  signal.  The  fluctuations 
in  load  are  on  the  order  of  magnitude  of  the  noise  level  present  in  the  signal 
and  it  is  not  possible  to  separate  the  two  effects.  An  additional  factor  to 
consider  is  the  system  compliance.  The  spring-in-series  used  in  these 
experiments  exhibited  some  hysteresis  effects  as  as  shown  in  Figure  8  which  is 
a  plot  of  spring  displacement  versus  applied  load.  The  spring  displacement 
and,  hence,  the  system  compliance  were  calculated  under  the  assumption 

6  TOTAL  =  6  SPRING  +  6  LOAD  LINE 

The  total  system  deflection  was  measured  with  an  LVDT  and  the  spring  deflection 
was  calculated  using  equation  (5).  In  all  calculations  of  Tappiie(j,  the 
stiffness  of  the  system,  KM,  was  assumed  to  remain  constant  and  its  value  was 
taken  as  the  slope  of  the  load-spring  displacement  curve  after  maximum  load 
since  this  corresponded  to  the  condition  of  the  system  when  the  instability 
occurred.  (Calculation  of  the  system  compliance  incrementally  over  the  course 
of  the  test  produced  large  oscillations  in  the  Tappne(j  curve  which  would  make  it 
difficult  to  make  any  conclusions  about  the  applicability  of  the  Ernst  analysis.) 
The  hysteresis  is  probably  due  to  friction  between  the  rollers  and  the  spring 
plates,  although  other  factors  may  have  contributed  as  well.  The  experimental 
procedure  assumed  that  only  one  parameter  would  be  varied  from  test  to  test, 
namely,  the  load  train  compliance.  Since  the  material  behavior  was  not 
consistent  from  specimen  to  specimen,  there  were  two  free  variables  interacting 
and  it  was  not  possible  to  conclusively  determine  the  effect  of  load  train 
compliance  on  the  amount  of  stable  crack  growth. 

One  additional  point  that  has  not  been  addressed  is  the  observation  of 
fracture  mode  transition  during  the  instability  event.  In  several  of  the 


tests,  the  fracture  mode  converted  from  ductile  tearing  to  cleavage  fracture 
after  several  hundred  mils  of  crack  extension.  Figure  9  is  a  photograph  of 
the  fracture  surface  of  specimen  GAF-5.  In  this  case  the  fracture  mode 
converted  to  cleavage  after  0.372  in.  (9.4  mm)  of  ductile  tearing.  This  is 
important  considering  the  fact  that  ductile  tearing  is  controlled  by  both  the 
stress  and  strain  in  the  vicinity  of  the  crack  tip,  while  cleavage  is  a  stress 
controlled  situation. 

CONCLUSIONS 

The  purpose  of  this  investigation  was  to  evaluate  tearing  instability 
theory  using  experimental  techniques  which  allow  constant  crack  length 
estimation  and  therefore  an  accurate  calculation  of  the  applied  and  material 
tearing  moduli  at  the  point  of  fracture  instability.  The  applied  tearing 
modulus  was  varied  from  test  to  test  using  a  spring  in  series  with  the  load 
train  of  the  testing  machine  to  attempt  to  produce  ductile  instabilities  at 
different  levels  of  crack  growth. 

Two  conclusions  can  be  drawn  from  the  results  of  this  investigation.  In 
most  cases,  fracture  instability  occurred  when  the  difference  between  the 
Ernst  ^appHed  formulation  and  Tmater^a2  was  on  the  order  of  10%.  The  effect 
of  changing  the  amount  of  stable  crack  growth  prior  to  fracture  instability  by 
changing  the  system  compliance  was  obscured  due  to  the  sensitivity  of  Tapp^ie(l 
in  the  vicinity  of  maximum  load  and  material  variability  in  the  A106  steel. 

It  was  also  observed  that  the  A106  steel  specimens  exhibited  large 
differences  in  the  resistance  curve  behavior  and  that  the  fracure  mode 
transitioned  from  ductile  tearing  to  cleavage  fracture  during  the  instability 
event  in  two  of  the  tests. 
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